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We study the lepton flavor violating (LFV) decays Z — l^l-^ — ''') ^^'^ 

framework of the minimal 331 model. The main contributions arise at the one-loop level 
via a doubly charged bilepton with general LFV couplings. We obtain an estimate for the 
corresponding branching ratios by using the bounds on the LFV couplings of the doubly 
charged bilepton from the current experimental limits on the decays li — >■ Ij^ and 
l~ — ^ ^j^^k^k' ^ bound on the bilepton mass is also obtained through the current 
limit on the anomalous magnetic moment of the muon. It is found that the bilepton 
contributions to LFV Z decays are not expected to be at the reach of experimental 
detection. In particular, the branching ratio for the Z — > /i^r^ decay is below the 10"^*^ 
level for a bilepton mass of the order of 500 GeV. 

PACS number(s): 13.38Dg, 13.35.-r 

1. Introduction 

In the standard model (SM), neutrinos are considered massless and thus lepton 
flavor violation (LFV) is forbidden at any order of perturbation theory. Even if the 
theory is extended with massive neutrinos, LFV transitions such as ii — ^ £j"f would 
be induced up to the one-loop level and would be strongly suppressed due to a 
GIM-like mechanism: it was found that BR(^ 67) ~ lO'^^ - IQ-''^ in the SM 
extended with non-diagonal lepton flavor couplings and massive neutrinos with a 
mass rrii, of a few eVs.^^^ Any signal of LFV would thus be a hint of new physics. 
However, recent evidences of neutrino oscillations and thereby a nonzero neutrino 
mass clearly point to LFV and have thus triggered the interest on the study of 
LFV decays such as £^ i^j, l~ ^pk^k^ and Z ifif {(ij = e, fi, r). 
Several theoretical extensions of the SM do predict such LFV transitions with a 
non-negligible rate. Currently there are stringent experimental constraints on LFV 
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muon decaysH BR(^ -> 67) < 2.4 x 10'^^^ BR(^ ^ 3e) < 1.0 x lO-^^J^ and 
BR(/LtTi eTi) < 3.6 x 10"".^ Even more, the bound on the /Lt — > ej rate is 
expected to be improved up the level of 10~^^ by the MEG experiment.^^ On the 
other hand, the current bounds on LEV transitions involving the r lepton are less 
stringent: BR(r ^ ej) < 3.3 x 10-^'£l BR(t -^ji-f) < 4.4 x 10-^,^1 BR(t 3e) < 
3.6 X 10-^1^ and BR(t e-e+^J.) < 3.7 x 10"^.^ As a matter of fact, the possibility 
that LEV transitions involving the last two lepton generations may be larger than 
those involving the electron has been widely conjectured in the literature. As far as 
LEV Z decay s are concerned, the most stringent experimental bounds were obtained 
at LEPM 

BR{Z-^e^H^) <1.7xl0-^\ (1) 
BR(Z ^ e^r^) < 9.8 X 10"^ (2) 
BR (Z ^ /i=Fr±) < 1.2 X 10"^ (3) 

The study of these LEV Z decays has been the source of great interest as 
the future international linear collider with its Giga-Z option would allow a 
yearly production of 10^ Z bosons,^^ which would open up the possibility 
for detecting some Z boson rare decays. Several predictions for the Z — > 
ifij decay width have been obtained in the framework of various SM exten- 
sion sM'4m^ ^^^^ such as supersymmetric theo- 

ries,13I2II2Slthe two-Higgs doublet model^^^the Zee-modelj^^lthe scalar triplet 
model,II21 the left-right sym metric 

model,!!^ 

top-color assisted technicolor,'^ the 
SM with massive n eutrino s, EHH effective theories, etc. 

The 331 model HHEH is an appeahng SM extension based on the SU{3)l x 
U{l)x gauge group. This model predicts new physics at the TeV scale and it is also 
attractive due to its peculiar mechanism of anomaly cancellation, which requires 
that the fermion family number is a multiple of the quark color number, thereby 
suggesting a solution to the flavor problem. A remarkable feature of the 331 model is 
the prediction of new exotic particles with masses bounded from above at the TeV 
scale due to theoretical constraints. Therefore, such a model could be confirmed 
or ruled out in a near future. Among the new particles predicted by the model, 
there are singly and doubly charged scalar bosons, exotic quarks of electric charges 
— 4/3e and 5/3e, singly and doubly charged gauge bosons, and an extra neutral 
gauge boson. The new charged gauge bosons are known as bileptons as they carry 
two units of lepton number. LEV can be induced at the tree-level in the scalar and 
gauge sectors. In particular, we will focus on the possibility that the new doubly 
charged gauge bileptons can give rise to LEV at the tree-level, which in turn can 
induce LEV Z decays at the one-loop level. Our aim is to present such a calculation 



^AU the experimental limits used in this work correspond to the 90% C.L. limits unless stated 
otherwise. 
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and obtain an estimate for the Z — )■ ^^^J branching ratios. To constrain the LFV 
bilepton couphngs we wiU use the current experimental bounds on the LFV decays 
ii — >■ and l~ — > i^i^^k ■ "^^^^ anomalous magnetic moment of the muon will be 
used to constrain the bilepton mass. 

The rest of our presentation is organized as follows. In Sec. [2] we present an 
overview of the minimal 331 model and consider the possibility of LFV mediated 
by the doubly charged vector bilepton. Section. [3] is devoted to the presentation of 
our calculation, whereas the numerical results and the analysis are presented in Sec. 
m The conclusions and outlook are presented in Sec. [3 



2. The minimal 331 model 

The minimal 331 model is based on the SUc{3) x SU{3)l x U{1)x gauge group. In 
this model , neutr inos are massless and the leptons are accommodated in antitriplets 
of 5C/(3)lP1^ 



z/i :(1,3*,0), 



(4) 



where i = 1, 2, 3 is the generation index and e^* is the complex conjugate field of e]^. 
Anomaly cancellation requires that the first two quark generations are represented 
by triplets of SU{3)l, while the third one appears as an antitriplet of SU{3)l. 
These multiplets are completed by three new exotic quarks (I?, S and T) with 
electric charges Qd,s = — 4/3e and Qt = 5/3e. For the purpose of this work it is 
not necessary a further discussion on the quark sector. 

The most economic scalar sector of the minimal 331 model requires three 
scalar triplets and one sextet of SU (3) One scalar triplet is necessary to break 
SU{3)l X U{l)x down to the electroweak gauge group, whereas electroweak sym- 
metry breaking (EWSB) requires the two remaining scalar triplets and the sextet. 
The latter is necessary to provide realistic masses for the leptons. This minimal 
Higgs sector has the following quantum numbers 



•/-y^^^o): (1.3,1); ^i=[l-)-- (1.3,0); <^2=(^^*'_j: (1,3,-1), 

(5) 

where $i = {(f>f,(f>i), with <i>j = it^^* for i = 1, 2, 3; $y = ("I>y+,"I>y) contains 
the would-be Goldstone bosons associated with the new doubly charged, , and 
singly charged, , bileptons; the real and imaginary parts of (fP correspond to one 
physical Higgs boson and the would-be Goldstone boson associated with the extra 
neutral gauge boson, Z' . In addition the scalar sextet is given by 
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H=\ h (1-6,0), (6) 

where T is a SU{2)l triplet 



~ \T+/V2 TO J ' 

whereas A~, p , and rj are singlets of SU{2)l with hypercharge —2, —4, and 
+4, respectively. 

The covariant derivative in the fundamental representation of SU{3)l x U{1)x 
can be written as 



D^ = d^-tg^W;-igxX^X^, (a = 1...8), (8) 

with A'' the GcU-man matrices and A'^ = \/2/3 diag(l,l,l). The first stage of 
spontaneous symmetry breaking (SSB) is triggered by the vacuum expectation value 
(VEV) of (pY, which breaks the SU{S)l x U{1)x gauge group down to SU{2)l x 
U{1)y- In this stage of SSB, the exotic quarks and the new gauge bosons acquire 
their masses. The bileptons appear in a SU{2)l x U{1)y doublet with hypercharge 
3 and are mass degenerate. They are defined in terms of the gauge eigenstates as 
follows 



(9) 



The gauge fields and Xfi mix to produce the extra neutral gauge boson, Z' , 
along with a massless gauge boson, B^, which is associated with the U{1)y group. 
They are given by 

Z'^ = ceW^^-seX^, (10) 
B^, = seW^ + ceX^, (11) 

where se = sm6,ce = cos9 and tan^ = gx/{V2g). The coupling constant associ- 
ated with the hypercharge group is g' = g sg / a/3. The remaining fields associated 
with the unbroken generators of SU (3) l are the gauge bosons of the SU (2) l group, 
which are denoted by for t = 1, 2, 3. 

After the first stage of SSB wc arc left with the SM with its particle content plus 
the new gauge bosons and exotic quarks, together with several scalar multiplets 
of SU{2)l: three doublets $i (i = 1, 2, 3), one triplet H, and various singlets. 
Electroweak symmetry breaking (EWSB) proceeds at the Fermi scale via the VEV 
of the SU{2)l doublets < >o= Vij\f2 {i = 1, 2). By simplicity it can be assumed 
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that < H >o= 0. In this stage, the SM particles acquire their masses and the 
bileptons and the Z' boson receive additional mass contributions. The extra mass 
terms for the bileptons, which arise from the Higgs kinetic-energy sector, violate 
the custodial SU{2) symmetry. Therefore, the bilepton masses split: 




This mass splitting is bounded by the hierarchy of the SSB: |TOy+ — mY++ \ < 3m^. 
On the other hand, by matching the gauge coupling constants at the first stage of 
SSB, it is found that 

which means that s'^{mz') has to be smaller than 1/4. It was found tha t this 
condition implies that the new Z' boson cannot be heavier than 3.1 TeV.l^l^^ 
From this result and the symmetry-breaking hierarchy u ^ vi,V2,V3, it is inferred 
that the bileptons masses are smaller than mz' /2 ~ 1500 GeV. 



2.1. LFV in the 331 model 

The possibility of LFV in the 331 model was first analyzed in Ref. |30] and more 
recently in Ref. |33] in a more general context. In the scalar sector, LFV can be 
mediated by the neutral and charged scalar bosons. However, these interactions are 
expected to be very suppressed due to the smallness of the Yukawa couplings. We 
will thus not consider LFV mediated by scalar bosons in our calculation. 

As far as the gauge sector is concerned, the neutral Z' gauge boson cannot 
mediate LFV at the tree-level: it turns out that the Z' boson couplings to the 
leptons are flavor universal, so this gauge boson cannot mediate LFV at the tree 
level as the rotation of flavor states to mass eigenstates yields a diagonal coupling 
matrix. Moreover, it is interesting to note that the Z' gauge boson has a leptophobic 
nature as its couplings to a lepton pair are suppressed by the \/l — 4s^ factor. 
On the other hand, the interactions between the bileptons and the leptons can be 
written in flavor space as 



£ = + -j^^'^n'^LY^ + H.c. (14) 

After a rotation to the physical states is performed {I'j^ ^ = Ul,r^l,r and v'j^ = 
Ulvl) we are left with 



C = + ^t^y^v''vLY+ + H.c. (15) 

where we introduced the unitary flavor mixing matrix = U]j Ul- The singly 
charged bilepton effects on the Z ^t^J decay will vanish due to the zero mass of 
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Fig. 1. Feynman diagrams for the Z — > i^f^J^ decay at the one-loop level in the 331 model. In 
the Feynman-t'Hooft gauge there is an additional set of Feynman diagrams obtained by replacing 
each bilepton gauge boson by its associated would-be Goldstone boson, but these diagrams can be 
neglected since their amplitudes are proportional to the lepton masses. 

neutrinos. We will thus only need to consider LEV mediated by the doubly charged 
bilepton. 

3. Analysis of the Z — >■ if^if decay 

To calculate the Z — > ^f^J decay we will neglect the masses of the outgoing lep- 
tons but the internal lepton mass will be retained. The transition amplitude will 
be calculated using the Eeynman-t'Hooft gauge. In our approximation, the contri- 
butions to the Z — ^ ^t^J decay arise from the four Eeynman diagrams shown in 
Eig.[TJ Although there is also a set of Eeynman diagrams obtained by replacing each 
bilepton by its associated would-be Goldstone boson, Gy , the respective amplitudes 
contain terms proportional to products of the lepton masses due to the form of the 
Gy^dj coupling, so these contributions can be neglected from the calculation. Apart 
from the Eeynman rule for the coupling of the doubly charged bilepton to a lepton 
pair, which can be extracted from Eq. (jlSp . we only need the Eeynman r ule for the 
vertex. The latter is given in the Eeynman-t'Hooft gauge by!^^ 

Z^{k)Y--{ki)Y++{k2) = -^gzYY({k2 - + (fc - fca - fci)^5„^ 

- (fc - fcl - k2)uga^}j , (16) 

where all the particles are incoming and gzYY = 1 ^ 4s^, with sw = sinOw and 
cw = cos 9w- 
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To obtain the transition amphtude for each Feynman diagram we used the 
method of Ref. |36l which is meant for processes involving vertices with complex 
conjugate fields, such as our lepton number violating vertices. The Feynman pa- 
rameters technique was used to solve the loop integrals. In the massless outgoing 
lepton limit, the Z{p) — > if {Pi)iJ (pj) decay amplitude can be written as 

iM {Z ^ £fef) = Fl'u{p,)j^PLv{p,)e"ip), (17) 

with P]^ = (1 — 7^)/2. The Ff^ function depends on the internal lepton mass and 
the bilepton mass. It can be expressed in the form: 



■,3 



3 



Fl'{xk,XY) ^ Tj^Y.^lVl,*I{xu,Xy), (18) 
I TT CW 



with 



I{X],,Xy) ^^fniXk^Xy). (19) 



We have introduced the notation xt — m^/m?^ and xy = "T-y/™^, with rrik the 
internal lepton mass and my the bilepton mass. The /„ function stands for the 
contribution of the nth Feynman diagram of Fig. [T] After some lengthy algebra, we 
obtain: 



fi{xk,xy) = 2gzyy {Syk + 2) (By - Bky) - {Sy^ + 2a;y - Xk) Cyu 

(20) 



'-By 

2 



f2{xk,Xy) = 2gR 



1 + {Syk + 2) [Bk - Bky) - \Bk + (1 + SykfCky 



2 

+ 2gLXkCky, (21) 
(/a + h){xk.xy) = ffL (1 - Bky - dBky) , (22) 

where gL,R = <?y i g^, with gy = —1/2 + 2s^ and 17^1 = —1/2 the vector and axial 
Z-lepton couplings. Note that the amplitudes of the two bubble diagrams must be 
combined to obtain the right limit for massless outgoing leptons. In addition, the 
following functions (the subscript denotes the corresponding dependence) have been 
introduced: 
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By = 2{1- TY arccotrr - log(a;r)) + A, (23) 



Bh = 2{l-Xk arccothAfe - log(xfc)) + A, (24) 
BkY = 1 + 7^ {xk log(xfc) - Xy log(a;y)) + A, (25) 

OYk 



dBkY = — — 3 yxY -^k- "^^kXY log J j ' (26) 

^'^^ = {s^^) + P [j^) - - ^(^-)' (27) 

Cvk = G{xk) — G{xy) — 4arccot {ty) arctan 



TY 



26 



Yk 



- log (j^f^r) i^y) + 2 Re [H{t+)] , (28) 

where 6Yk = xy - Xk, A± = ^ (1 ± A^), Afc = a/I - 4a;fe, r+ = ^ (1 + iry) and 
Ty = ^/4x^^^. A stands for the usual ultraviolet singularity in dimensional regu- 
larization. Furthermore 

F(A) = log(A) log ( - log(A - 1) log (^^] + Li/ ^ 



C„ = log (S|±i^)+L«(^). (30) 

^W = "<7Ti^)-"=(7^)- 

Although each Feynman diagram is ultraviolet divergent by itself, all the divergences 
cancel each other out. This becomes evident when we write qzyy = —2gV = 
-{9l+9r)- 

We would like to point out that to cross-check our calculation, wc made an 
alternative evaluation via the Passarino-Veltman method, using the unitary gauge 
and without any approximation. The resulting amplitudes are rather lengthy to 
be included here, but numerical evaluation showed a nice agreement between our 
approximate result and the exact calculation. 

The Z £fej decay width is given in the massless outgoing lepton limit by 

nz ^ ifif) = ^ {\Fl^\' + \Fir) . (32) 
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Fig. 2. Feynman diagrams for the £^ — t^jt^j. i^. decay at the one-loop level in the 331 model. 



To obtain an estimate for the Z — > if£j branching ratios, we need to analyze the 
bounds on the mixing matrix and the bilepton mass. Below we wih examine 
the bounds obtained from the £j ^t^k^k ^^'^ ~^ ^i'^ decays together with the 
muon anomalous magnetic moment. 



4. Numerical analysis 

4.1. Bounds on the bilepton LFV couplings 

The three-body decay £~ ^^^k^k proceeds at the tree-level through the Feynman 
diagram of Fig. [2l Its decay width can be obtained straightforwardly in the limit of 
massless outgoing leptons. It is given by: 



4- '^ 

9 m° 



3 2ii7r3mt 



■w^ki'mi' + w,";]') 



(33) 



The experimental 90% C.L. limits on the LFV decays /i — e+e e pl r 



e + e — e— j^and r — > /i+e e^^translate into the following bounds with 90% C.L. 



\vK\H\v,l\' + \vy\') 



ml 



\v,\\H\v,l\' + \v,\\') 



< 2.57 X 10"^" GeV" 



< 5.14 X 10"^^ GeV~'', 



InYl' (1^231' + 1^321') 



ml 



(34) 
(35) 
(36) 



We used the fi and t mean lifetimes given in Ref. [37l Following Ref. [30j we will 

parametrize the mixing LFV matrix as V^j ~ Sij + ay C^^'J , where = —aji and 
6ij = —9ji are mixing angles and CP-violating phases, respectively. Therefore we 
obtain the foUowint limits with 90% C.L.: 
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Fig. 3. Feynman diagrams for the ii — > ij^ decay at the one-loop level in the 331 model. 



< 1.28 X 10-20 GeV-^ (37) 



< 2.57 X 10-^2 GeV"'', (38) 

my 

< 2.64 X 10-12 GeV-'*. (39) 

rriY 

As far as the — ?> tj^ decay, it proceeds via the triangle Feynman dia- 
grams shown in Fig. [31 The bubble diagrams does not contribute in the limit of 
massless outgoing lepton. Using the same scheme described above to obtain the 
Z — >■ tftj amplitude, we obtain the ii — >■ £^-7 amplitude, which can be written as 

M{i, ^ ^,7) = u{p,)<j^^ (r^ + fl^j') u{p,)e^{p), (40) 

nii + nij V ^ / 

where q = Pi — Pj and pi (pj) are the photon and the incoming (outgoing) lepton 
4-momenta, whereas 



/Ji. - ^ I r (;! + ^) i (vj vi- ± v-r V?.-) (41, 

with Ml ^x{l-y + (1 - x)^k , M2 = x{^k -yS,i) + {l-x) and ia - m^/mf,. We 
have dropped any terms independent of the internal lepton mass, which cancel due 
to the unitarity of the matrix. Furthermore, if one neglects the lepton masses 
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in the Mi and M2 coefficients, it follows that 



fiA = E ^^^^ ± ) (42) 



3 



267r2 



fe=l 



which coincides with the result previously found in Ref. 1301 In the massless outgoing 
lepton limit, the — )■ ij^ decay width can be written as 

nh^hi)-'-^{^m + \nf)- (43) 

Since the sum in Eq. (j42p depends on the internal lepton mass and the [i" — !■ e+e^e^ 
decay gives a strong constraint on we can neglect the electron and muon terms. 
It follows that the experimental 90% C.L. limits on the decays /i — > 67, ^ r — 67,^ 
and r — )■ /Lt^EI translate into the following bounds with 90% C.L. 

'^i3l'l^23l' < ^ 29 X 10-21 GeV-4, (44) 



my 



< 2.35 X lO-'-" GeV-^, (45) 

< 3.14 X 10-" GeV--*. (46) 



As far as the lepton anomalous magnetic moment ag — {g~ 2)/2 is concerned, it 
can receive new contributions in the 331 model from the scalar and gauge sectors. 
The doubly charged bilepton contributes through triangle diagrams analogue to 
those of Fig. [H Its contribution is given by 

2 3 

" ■ - ^ ^ |V,IP (9mfe cos{2e,u) + (47) 



^ k=\ 

If the flavor mixing matrix is diagonal, ai becomes 



(48) 



which agrees with the result obtained before in Ref. [38l We also observe that the 
contribution to from the scalar bosons as we ll as the singly charged bilepton and 



the extra neutral gauge boson are subdominantP^, so the 331 model contribution, 
a^^i, can be assumed to arise mainly from the doubly charged bilepton. The current 
experimental limit on is not useful to bound the LEV couplings but it can 
constrain the bilepton mass my. For the theoretical and experimental values of a^, 
we will use the most recent data quoted in Ref. |37l The discrepancy between the 

theoretical SM contribution, a^^, and the world average, a^"^' , of experimental 

39 • ■ 1 

measurements -"^ is given by: 
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Aa^ = a^^'P- - af^ = 255 (63) (49) x 10~", (49) 

where the e~e+ — )■ 7r~7r+ data from BABAR were used to evaluate the hadronic 
contribution to a^'^ Although this discrepancy is about 3.2 standard deviations, 
it is not yet conclusive since there is stih a considerable discrepancy between the 
various evaluations of the hadronic contribution. If Aa^ is ascribed to the doubly 
charged bilepton, we get the following bound on my with 95 % C.L.: 

my > 421 GeV. (50) 

Let us assess how this bound compares with other indirect bounds. The very strin- 
gent bound my > 800 GeV was obtained from muonium-antimuonium conversion. 
I^This bound, which would rule out the minimal 331 model, is based on the assump- 
tions that the matrix is flavor diagonal and the scalar sector of the model does 
not contribute significantly to muonium-antimuonium conversion. Another strin- 
gent bound, TOy > 750 GeV, arises from fermion pair production and lepton-flavor 
violating processes. ^21 These bounds can be evaded if one considers an extended 
Higgs sector or less restrictive assumptions.'^ We will rather consider a bilepton 
mass of a few hundreds of GeV to obtain an estimate of the Z^ifiJ branching 
ratios. 



4.2. The Z — >■ i^if branching ratios 

Due to the unitarity of V"^ , i.e. J^k ^ik^kj" — ^ij^ ^^^"^ neglecting imaginary phases, 
we can write 



BR(Z -> ijtf) = A V^^VJ^ {I{x2,xy) - I{xi,xy)) 

+ VlVl (/(X3,xy) - /(xi,xy)) |', (51) 

with A = YiTxTz (2^1^17) ■ Notice that according to Eq. ([37)) . is strongly 

constrained, so one can write 

BR(Z ^ eT^±) ~ A l^i^l' |t/2^|' |/(x3,Xy) - /i(x3,Xy)|2 , (52) 

BR(Z-^eTT±) ~ A|Fi^f |/(x3,a:y)-/(xi,a:y)p, (53) 

BR(Z^^Tt±) ~ A|F2^f |/(a;3,a;y)-/(a;2,a;y)|'. (54) 

Numerical evaluation together with the bounds (j44 |) - (|46l) give the upper bounds 
on the LEV Z decays shown in Table [1] for my = 100 GeV and my — 500 GeV. 
While the loop amplitude magnitude decreases for larger values of my , the bounds 
on the LEV matrix elements loosen up. This explains the fact that the bounds on 
the LEV Z decays are slightly weaker for larger my. However, the bounds obtained 
for my — 1000 GeV are of similar order of magnitude than those obtained for 
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Table 1. Upper limit with 90%C.L. on the bilepton con- 
tribution to LFV Z decays. The bounds given in Eqs. 
II44I I-I I46I I were employed. 



BR 


my = 100 GeV 


my = 500 GeV 


BR(Z -5- eT/i±) 


3.2 X 10"^" 


l.lx-W 


BR{Z -5> eTT±) 


6.37 X 10-" 


4.57 x 


BR(Z ^ M^T±) 


6.56 X 10"" 


2.26 X 10"^" 



my = 500 GeV. Our results indicate that the bilepton mediated Z decays would 
be far from the reach of detection, though the corresponding branching ratios are 
of the same order of magnitude than in other SM extensions. For instance, in the 
framework of the Zee-model it was found that BR(Z e^^^) < 4.2 x 10~^^, 
BR(Z e=FT±) < 1.1 X 10-8, and BR(Z ^ /x^r^) < 1.6 x lO-^^ for typical values 
of the model parametersfl^ also, in the SM enlarged with massive neutrinos with 
masses of the order of a few dozens of GeV, the upper bound BR(Z fi^T^) < 
10"" was obtained.lSl 

5. Final remarks 

We have investigated the Z — > iftj decay in the framework of the minimal 331 
model. We focused on the contributions mediated by the doubly charged bilepton 
and estimated the bounds on the LFV bilepton couplings from the experimental 
constraints on the decays £i — >■ and £~ — >■ £j'£k^k ■ results indicate that 
these contributions seem to be far from experimental detection. The smallness of the 
decay rates can be explained mainly from the fact that the LFV bilepton couplings 
are strongly constrained by current experimental data. 
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